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INTRODUCTION 



A. THE PACKET SWITCHED NETWORK CONCEPT 

A new technique for data communications that has evolved 
over 10 years is called PACKET SWITCHING. In general, commu- 
nication networks may be conveniently divided into three 
types: CIRCUIT SWITCHING, MESSAGE SWITCHING and PACKET 

SWITCHING. Both message and packet switching uses a tech- 
nique known as store and forward transmission. 

A circuit switching network provides service by setting 
up a total path of connected lines from the source to the 
destination of the call. This complete circuit is set up by 
a special signaling message that threads its way through the 
network, seizing channels in the path as it proceeds. After 
the path is established, a return signal informs the source 
that data transmission may proceed, and all channels in the 
path are then used simultaneously. 

The entire path remains allocated to the transmission, 
and only when the source release the circuit will all these 
channels be returned to the available pool for use in other 
paths. Circuit switching is the common method for telephone 
systems . 

In message switching, only one channel is used at a time 
for a given transmission. The message first travels from 
its source node to the next node in its path, and when the 
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entire message is received at this node, then the next step 
in its journey is selected. If this selected channel (link) 
is busy, the message waits in a queue, and finally when the 
channel becomes free, transmission begins. Thus, the message 
"hops" from node to node through the network using only one 
channel at a time, possibly queueing at busy channels, as it 
is successively stored and forwarded through the network. 

Packet switching is basically the same as message switch- 
ing except that the messages are decomposed into small equal 
pieces called packets, each of which has a minimum length. 

These packets are numbered and addressed and make their way 
through the net independently of each other. Thus, many 
packets of the same message may be in tranmission simulta- 
neously, giving one of the main advantages of packet switching. 

With packet switching systems, information is exchanged 
in the form of short packets. A packet-switched network can 
handle several different types of traffic concurrently. These 
include HIGH-THROUGHPUT traffic, for example, the transmission 
of large data files between computers, for which accuracy and 
high average data speed are the most important performance 
requirements; LOW-DELAY traffic, for example, interactive 
communication between a person at a terminal and a remote 
computer, for which accuracy and low average message delay are 
important; and REAL-TIME traffic, for example, packetized 
speech for which the performance of circuit-switched connec- 
tion must be approached by maintaining a relatively constant 
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data speed, but for which extreme accuracy is not important 
owing to the redundancy of the information. 

The packet switched network is designed primarily for 
computer to computer communication. It has a much more 
rapid response which matches the internal behavior of compu- 
ters and handles information in much the same way as does a 
computer. At the same time it can readily match the speed 
of attached computers to that of the terminal users, by virtue 
of its internal storage. 

The prime purpose of store and forward packet switching 
is to enable communications resources to be used effectively 
and in such a way that they may be shared by many users opera- 
ting in an intermittent fashion, giving each user a rapid 
response from the communication network just at the instant 
when thxs is required. 

If there is need for transmitting a long continuous stream 
of data, then a circuit switched connection makes good sense. 
On the other hand, if the data flow is bursty, then some form 
of resource sharing can be used to great advantage; packet 
switching is an effective choice here. 

Since packets are stored as they pass through switching 
nodes, it is possible to conduct speed, format and code con- 
version during the switching process. Another feature of 
packet switching is its ability to adaptively select good 
paths for packet journeys as a function of the network 
congestion. 
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Besides providing small network delays, packet switching 
has the desirable feature of rapidly handling small messages 
in spite of the presence of long messages that may be in 
transport at the same time, this is because of the decomposi- 
tion of long messages into packets. Another useful property 
of this decomposition is that the nodal storage requirement 
is reduced. 

In evaluating packet switching networks, the delay, 
throughput, cost and reliability are important measures. 

Theoretical studies have been directed to the queueing and 
network flow problems in general and more specifically to such 
problems as delay analysis, route assignment, topological 
design and flow control, etc. 

The topic chosen in this study is centered around optimal 
flow and minimum capacity assignment and delay analysis for 
packet-switched networks. The networks under consideration 
have a relatively complicated structure with a large number 
of source-destination node pairs. 

3, THE ROUTING PROBLEM 

A message routing procedure is merely a decision rule 
that determines the node a message will next visit in its 
path through the network. The objective of the routing 
procedures is to transport packets on a minimum delay path 
from source to destination. 

In discussing routing policies for networks, an important 
distinction must be made between static and adaptive policies. 
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This distinction depends on the environment in which a policy 
is designed to operate. If network topology is not subject 
to changes (due to failure, modifications, growth) and traffic 
inputs are stationary, then the optimal routing solution is a 
static solution consisting of a set of fixed paths between all 
node pairs. The traffic between each source and destination 
pair may be distributed on several paths simultaneously in 
well defined proportions, where the proportions are fixed in 
time . 

In a real network environment, however, the topology 
changes with time and user traffic requirements tend to fluc- 
tuate more or less rapidly. To minimize delay it is then 
necessary to implement an adaptive routing policy that can 
react and adjust to various changes. The adaptivity of a 
policy can be measured in terms of its response time. A 
reasonable procedure is to use a periodically refreshed 
static routing solution. Indeed, there is a continuum of 
solutions between the two extremes, static and adaptive, 
characterized by different response times and used for 
different applications. 

Beside their use in operational networks, static routing 
policies find an important application in the network design 
process. During this process, for a given traffic pattern a 
prediction of the throughput and delay performance of a 
given topolocy is needed. The routing policy clearly has a 
major impact on such performance. Most routing implementa- 
tions are adaptive, and unfortunately the analysis of adaptive 
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routing policies is an extremely difficult task. To simplify 
the problem, the traffic pattern is usually approximated with 
a stationary pattern, and the routing policy with a static 
routing policy. 

1 . Static Routing 

a. Routing Table Representation 

A static routing policy is represented by a set 
of routing tables, one for each node, indicating how the 
traffic arriving at that node must be routed on the outgoing 
links, depending on its final destination. The routing table 
for node i is an Nx L matrix Pi{n,k) , where N is the number 
of nodes in the network and Pi{n,k) is the fraction of traffic 
directed to node k which, upon arrival at node i, is routed 
through neighbor node n (Fig. 1.1). By the definition. 

Pi (j,k) = 1 

Where n is the number of neighbors of node i. 

The actual distribution of the incoming traffic 
to outgoing links may be done randomly using as probability 
weights the values Pi(j,k). If for any node i there is only 
one permissible outgoing link from i to any destination k, then 
Pi(j,k) = 1 for neighbor node j. 

0 otherwise 

Such a routing policy is referred to as a single path 
routing policy, since only one path is used from any node i 
to any destination node k. In general, the optimal static 
routing solution is a multipath solution, allowing for the 
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simultaneous use of several routes in order to minimize 
delays. The routing table representation may also be used 
for adaptive policies. 





Figure 1.1. Routing Table. 

b. Problem Formulation 

In packet switched network, messages are seg- 
mented into small packets, and each packet travels from its 
source to its destination via a set of intermediate nodes. 
While awaiting transmission to the next node in the route, 
the packet must be stored until the link to that node is 
free. Thus, at each node there are several queues, one for 
each output channel. 
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Packet flow requirements between nodes arise at 



random times; therefore, link flows, queue length, and packet 
delay are random variables. A static routing problem can be 
defined as the problem of finding the static routing policy 
which optimizes the flow and minimizes the average time delay. 

Considering the relationship between routing policy 
and link flows, the problem can be formulated for solution of 
optimal flow and minimum capacity with respect to the given 
requirement as follows, 
given : (a) topology 

(b) requirement matrix R 
minimize: alpha 

subject to: (a) conservation of flow 

(b) flow £ capacity x alpha 
Where alpha is a nonnegative scale factor. 

The mathematical model and solution techniques are 
discussed in detail later. 

2 . Adaptive Routing 

The adaptive routing problem consists of defining a 
procedure that dynamically updates routing tables according 
to changes in the network. This procedure will include the 
acquisition of the status of all neighboring nodes, which a 
packet associated with the source and destination node pair 
could visit next after being processed at the current node. 

The neighboring nodes are those with direct channels linking 
to a given node . 
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The structure of the routing table and the ways the 
information contained therein are to be used vary according 
to the individual routing procedure. For example, if single 
path routing is used, only one path links each source and 
destination node pair; therefore, only a single neighboring 
node associated with each node pair is listed in the routing 
table . 

The design of routing is a process of determining the 
structure of routing tables and specifying the procedures of 
using them. The choice of routing depends on many factors 
such as network topology, available facilities, throughput 
and delay requirements, and so on. Adaptive routing tech- 
niques are ruled out in this study. 
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II. OPTII-IAL FLOW AND MINIMUM CAPACITY ASSIGNMENT 



The optimal flow is the best selection of paths between 
source and destination, given the traffic requirements and 
the network configuration. The definition of best paths may 
vary depending on the nature of the traffic. In this study, 
the best paths will be regarded as the paths of minimum 
"distance", where distance can be interpreted in various ways, 
for instance as time delay along the links. 

To investigate average time delay of the networks, two 
different approaches are discussed using classical optimiza- 
tion procedures . 

A. NETWORK MODELING 

In a distributed packet-switched network, packets are 
transmitted over a collection of switching nodes in tandem. 

In modeling a packet switched network, the objective is 
usually to optimize the total flow of the networks, to 
minimize the capacity of each link, or to maximize the 
utilization of the network. The problem may be formulated 
as follows. 

given topology 

requirement matrix R 

capacity of each link (arbitrary) 

minimize; alpha 

subject to: (1) conservation constraints. 
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IfI (k) 



Ri (k) 



lF£(k) 

outgoing incoming 

links links 

from node i to node i 

(2) capacity constraints 
Ft (k) £ alpha x Ct 

(3) positive constraints 
iFt (k) ^ 0 



input 
at node i 
for node k 



Where F£(k) is the flow on link destined for node k. We 
next look at each constraint further in detail. 

1 . Conservation Constraints 

At each node i, the total outgoing packet flows 
destined for node k equal the incoming flows plus the input 
to be transmitted to node k. 



iFe(k) 

outgoing 

links 

from node i 



iF^(k) 
incoming 
links 
to node i 



Ri(k) 
input 
at node i 
for node k 



where 

k = destination node (k = 1,2,3...N). 

N = number of nodes. 

L = number of links. 

i = source node (i = k, i = 1,2, 3. ..N). 

(! = link between two nodes {t= 1,2,3...L). 
F£(k) = flow on link destined for node k. 
Ri(k) = input destined for node k at node i. 
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2 . 



Capacity Constraints 



The total flow destined for node k on a certain link 
should be less than or equal to the capacity of the link i. 
Since linear programming algorithms have difficulties dealing 
with inequalities, slack variables are added to the constrains 
to absorb unused resources and thus to force equalities to 
appear . 

When slack variables are introduced, the problem in 
standard form becomes. 

F£(k) - alpha X c£. + = 0 

k=l 

where 

Ct = capacity of the link. 

alpha = constant to optimally minimize the link 
capacity and total flow on link. 

S£ = slack variables. 

3 . Positive Constraints 

If we consider the network to have unidirectional 
links, the bidirectional flows between two nodes are divided 
into two one-way links in opposite direction, and the flow on 
each link is positive. 

F£ (k) > 0 

Using matrix notation, the model simply may be 
expressed as a linear programming problem. 

minimize; alpha 

subject to: AX = R 

X >0 
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B. MODEL PROGRAxMMING 

In order to use a readily available Mathematical Program- 
ming System such as MPSIII, which can solve linear programming 
problems with up to 4000 rows and theoretically with an un- 
limited number of variables, the major programming problem 
with large scale networks is how to generate the input data 
for the MPSIII program. In other words, how do we generate 
the matrix form of the model. 

1 . Node to Link Incidence Matrix 

For any network processing to take place in the compu- 
ter, the network structure must be represented in some machine 
understandable form. A wide spectrum of different represen- 
tation methods are available. Most of these are based on an 
incidence matrix scheme. A matrix form is convenient for the 
analysis of networks since the resulting matrices are in a 
format suitable for mathematical analysis. 

The network of Fig. 2.1 will be used to illustrate the 
method. The node to link incidence matrix E(i,£) describes 
the links connected to the node such that 

E(i,£) =f 1 if link t is outgoing from node i. 

-1 if link Z is incoming to node i. 

^ 0 if link Z is not connected to node i. 

The incidence matrix for input at node i destined for node k 
is a modification of the node to link incidence matrix E(i,£) 
according to the following rules: 
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Ek(i/ ) 



0 for columns corresponding to links 



= r 

outgoing from node k, 
delete the row if i = k 
^E(i,£) otherwise 
For the illustration network (Fig. 2.1). 
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E(i,£) is a 3x5 matrix and Ek(i,£) is a 2x5 matrix in size. 




Figure 2.1. Illustration Network. 

2 . Matrix Representation 

Matrix representation is a convenient tool for 
programming. The complete matrix representation for the 
illustration network is shown in Appendix A. 

The general model matrix AX = R may be represented 
as follows: 
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where 



Ek = (N-1) X L incident matrix. 

Ik = L X L identity matrix with zero elements 

corresponding to zero columns of the Ek matrix. 

I = L x L identity matrix. 

Ri(k) = requirement matrix which is the inputs to be 

sent from source node i to destination node k. 

The size of the model in matrix representation is; 

A = ((N-1) X N + L) X ( (N+1) X L + 1) 

X = ( (N+1) X L + 1) X 1 

R = ( (N-1) X N + L) X 1 

For example, the network consisting of 20 nodes and 
40 links has the matrix which is A = 420 x 841, X = 341 x 1 
and R = 420 x 1 in size. 

3 . Programs 

To generate the MPSIII input data and to solve the 
model, requires two programs which are listed at the end of 
this study (Appendix C and D) . 

The first one is the data (model) generation program 
(Appendix C) written in the Fortran language. The output of 
this program uses special notation to designate the flow of 
the link connecting node pairs, destination node of the flow, 
flow variables, slack variables and input at the node. These 
notations are composed of a letter which is one of C, L, X or 
S and 7 numerical digits. 
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The notation Ci, Lt, X£, and S£ imply conservation 
constraints, capacity constraints, flow variables, and slack 
variables respectively. The rest of the six numbers following 
Ct, Lt, Xt and St indicates link connecting node pair with 
first 4 digits and destination node or slack variable number 
with last 2 digits. For example, the notation C1000201 indi- 
cates the conservation constraint row which is the input at 
node 2 destined for node 1, L1020100 the link capacity row for 
the link connecting between node 2 and node 1, X1020101 the 
flow variable on the link connecting between node 2 and node 
1 destined for node 1 and SlOOOOOl the first slack variable. 

The second program is the MPSIII program developed by 
Management Science Systems, Inc. It provides elaborate con- 
trol language for the formulation of solution strategies for 
mathematical programming problems. 

Individual instructions in this language bring in 
quite elaborate sections of code designed to execute the step 
as efficiently as possible. 

Using the special notation discussed above, there is 
a limitation on the size of the problem which can be solved, 
because data names must consist of 8 characters maximum. 

This means that the maximum network size must be less than 
100 nodes, if the computer capacity allows. 

C. SOLUTION TECHNIQUES 

As mentioned before , the classicial optimization proce- 
dures to solve the investigation model provide solutions 
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which are the optimal paths and minimum link capacities for 
the network. Different optimal solutions may be expected 
depending on different methods of optimization. The two 
different approaches, each involving iteration of the linear 
programming solution procedure are discussed below, and illus- 
trated in terms of one example network problem. These 
approaches are based on the characteristics of linear pro- 
gramming problems. 

1. Successive Saturation Approach 

The standard procedure of the MPSIII program listed 
at the end of this study (Appendix D and E) produces an 
"optimal" solution for alpha from the output data of the data 
generation program in the first run. But the corresponding 
flow solution is not a good solution in terms of the average 
time delay in the packet-switched network, when it is ana- 
lyzed by the Kleinrock's delay analysis model (Ref. 4) . The 
problem is that only the flow thru the saturated links is 
optimum. 

In order to reduce the time delay. Iteration runs are 
attempted until all links are saturated. The detailed itera- 
tion method for the example network (Fig. 3.1) which has 5 
nodes and 6 links, each with original link capacities 10 (to 
simplify the problem) is discussed below. 
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Figure 3.1. Example Network. 

The ouput of MPSIII program is divided into row 
section and column section (Ref. 9 and Appendix D and E) . 

The objective function value (which is the same as the alpha 
value) and the links saturated with the alpha value in the 
capacity constraint rows indicated with the first letter 
notation L can be read, and the flows on each link and the 
alpha value in column section. In the example network, for 
the first iteration the saturation level (alpha value) is 
0.25, and the saturated links are L1010300 (link 3) and 
L1040300 (link 7) . Their activity and slack activity levels 
are all zero, but their dual activity level is not (pp 48) . 

For the next iteration, the changed input datas are 
that the original capacity multiplied by the saturation level 
of first run are moved from the alpha column to the right 
hand side (RHS) column. In the example, the new link capaci- 
ties 2.5 (0.25 X 10) are moved to RHS (pp. 50 ) , and the 
problem is solved again. Repeating this iteration procedures 
until all links are saturated and moved to RHS, the objective 
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fxinction value becomes zero. When the objective function 
value is zero, the final flow solution is the desired solution 
of this successive saturation approach. In the example case, 
the six iterations are required to solve the problem. These 
whole procedures except programs and input datas from 2nd run 
for the example network are listed at the end of this study 
(Appendix D) . 

2 . Max- Slack Approach 

To simplify the solution of the above successive satu- 
ration approach, another approach called the max-slack approach 
is considered. It requires only two iteration of the program, 
but is otherwise quite similar. Like in the successive satu- 
ration approach, after the first run, all saturated link 
capacities are multiplied by the alpha value of the first run 
are moved to the RHS , and in addition, all slack variable 
columns are summed as the objective function in the data card 
for the second run. 

In the example network, eight slack column variables 
are added on the objective rows, and all new capacities 2.5 
(0.25 X 10) are moved to RHS. This is clearly shown on the 
output of control language "picture", and the program, input 
data and the output are also listed at the end (Appendix E) . 

D. DELAY ANALYSIS 

The average time delay of a network is the average time 
a packet spends in the network traveling from its source to 
its destination nodes. To obtain a tractable expression for 
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average time delay T in terms of the capacity of link i, it 
is necessary to make simplifying assiamptions . Packet lengths 
are assumed to be exponentially distributed, and are re-chosen 
with statistical independence at each node to form a Poisson 
process . 

These assumptions ’ known as the 'Independence assumption', 
introduced by D. Kleinrock (Ref. 4) , make each queuing prob- 
lem independent. In the analysis, each queue in a packet 
network is assumed to be an M/M/1 system with arrival rate 
of gZ and to be independent of each other. 



The average delay for link t is given by the waiting time in 
the queue as 



T 



I 



1 

uC^-G 



t 



Where G£ is the packet traffic in the link and 1/u the 
average packet length in bits. To form a suitable average 
over all the queuing processes, the T are weighed by Gt/r 
where r is the total input packet rate to the network. In 
this way, the total delay suffered by all packets per second 

of network operation I Gt*T is divided by the total number 

i 

of packets carried by the network per second. 

The total average time delay per packet through the net- 
work is then given by 

T = L V , L y 

^ 1=1 ^Z=l 
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Where Gt/u = Ft, and L is the total number of links in 
the system and 1/u is equal to the average flow on link t. 

It should be noted that the value of Ct (the minimum required 
capacity) must be greater than the total of all the flows on 
link t. In other words, the value of Ct must be greater than 
the value which is the first run saturation level multiplied 
by the arbitrary capacity C for the first run. 

The value of T includes both of the waiting time and the 
service time. The waiting time is subject to the interference 
of all other traffic within the network that consists of the 
data traffic as well as the control traffic. This average 
time delay can be used to analyze the relative goodness of the 
routing strategies considered in this study. 

In deriving this basic formula, a number of factors (such 
as processing time and propagation delay) are neglected. In 
any realistic network, these variables as well as others must 
be considered in the analysis of the networks. According to 
the above formula, the average time delay depends only on the 
aggregate flow of the links. The reduction of the flow Ft of 
each link results in the reduction of the average time delay 
for the network which has constant input packet rate and link 
capacities . 

For the example network (Fig. 3.1), the arbitrary capa- 
city is 10 and the value of Ct must be greater than 2.5 
(0.25 X 10) . The average time delay for the capacity F£=10, 
and the input (requirement matrix) of 2 units from node 1 to 
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node 4 and 5 units from node 2 to node 3 are Tss = 2.12/r 
for the successive saturation case and Tms = 2.083/r for the 
max-slack case. To make these calculations easy, the data 
summarizing the final output run are listed in Table I. We 
see that there is not much advantage of one procedure over 
the other in terms of performance, whereas the max-slack 
solution requires only two iterations, rather than 6 for the 
successive saturation algorithm. 
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TABLE I 



LINK FLOW FOR THE EXA^IPLE NETWORK 



( ) ; Max-Slack 



LINK 




DESTINATION NODE 


TOTAL 


(Nodes) 


1 


2 3 4 5 


FLOW 


1-2 




(2) 


(2) 






2 


2 


2-1 




(2.5) 


(2.5) 






2.5 


2.5 


1-3 




(2.5) 


(2.5) 






2.5 


2.5 



3-4 



2-4 


(2.5) 

2.25 




(2.5) 

2.25 


2-5 


0.25 


(2) 


(2) 






2 


2.25 


4-3 


(2.5) 




(2.5) 




2.5 




2.5 


5-4 


0.25 


(2) 


(2) 






2 


2.25 


Req Mat 


O 
1 — \ 
It 

U 


R^(4)=2, R2(3)=5 
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III. EXPERIMENTAL RESULTS AND CONCLUSIONS 



A. EXPERIMENTS AND COMPARISON OF THE RESULTS 

The experiment was conducted for two experimental networks 
to analyze both approaches, in addition to the example network. 
These networks are specified in Appendix B. The networks were 
used in the experiment with requirement matrices of 15 units 
from node 1 destined for node 5, 9 units from node 7 destined 
for node 2 for the 9 nodes/36 links network (Fig. B.l) ; and 
15 units from node 2 destined for node 12, 4 units from node 
5 destined for node 6, and 10 units from node 10 destined for 
node 1 for the 13 nodes/60 links network (Fig. B.2). Both the 
successive saturation and max-slack approaches were used for 
the successive saturation approach, the 21 iterations are 
required in this experiment. 

The results and their comparison is listed in Tables II, 
III, and IV. These are obtained by writing each flow variable 
value from the final output run on each link of the networks, 
and summing these flow variable values to obtain the total 
flow on the link (Appendix F and G) . 

As shown in Tables II and III, the number of iterations 
required and the total link usage (utilization) are much 
greater for the successive saturation approach, but fewer 
large link capacities are required than in the max-slack 
approach. The average time delay depends on the capacity 



34 



RESULTS FOR THE 9/36 NETWORK 
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RESULTS FOR THE 13/60 NETWORK 
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value Ci (Table IV) . When the capacity value Ct is close 
to the minimum required capacity, the time delay is less for 
the successive saturation approach than for the max-slack 
approach in both networks, and vice versa. It is significant 
for the large scale network, and for large enough link capa- 
cities CZ , that the time delay difference becomes negligible. 

Also, in the process of this study, it was noted that both 
approaches yield the same aggregate flows in the case of very 
simple networks (nodes less than or equal to 4) , and for com- 
pletely symmetric networks and requirement matrices (input and 
capacity matrices) . For example, the time delays are the same 
for both approaches with 2 units input from the outmost nodes 
to the directly opposite nodes respectively (1 to 13, 2 to 12, 
5 to 9, 12 to 2 and 13 to 1) for the 13 node network, because 
all links total flow are 3/4. The detail experimental results 
are listed in Appendix H. 

B. CONCLUSIONS 

This is an elementary study of the routing design problem 
for a packet network: given a traffic requirement matrix, 

minimize the average time delay per packet, subject to finding 
a feasible flow for a network with fixed topology and link 
capacities . 

The general performance characteristics and advantages 
of the two approaches are investigated for small and simple 
networks. This illustrates one way to solve the design 
problem under limited conditions. 
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TABLE IV 



TIME DELAY FOR THE CAPACITIES 



CAPACITY 


4 


4.2 


4.5 


5 


10 


15 




SUCC 




244. 1/r 


109. 5/r 


59 .9/r 


10.74/r 


5.82/r 


9/36 


SAT 














Net 


MAX 

SLACK 




28.56/r 


116. 6/r 


59. 5/r 


10 . 41/r 


5.73/r 




SUCC 


169. 9/r 


115. 6/r 


75. 7/r 


60. 1/r 


15. 34/r 


8. 9 3/r 


13/60 


SAT 














Net 


MAX 

SLACK 


335. 7/r 


188. 3/r 


114. 2/r 


69. 3/r 


14.29/r 


7.99/r 
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APPENDIX A 

MATRIX REPRESENTATION FOR ILLUSTRATION NETWORK 



A 


X 




R 


0 10 0 1 




F^d) 




R2d) 


0 0 0 1 1 




F 2 O) 


= 


Rad) 


10 1-10 




F3(1) 




Rl(2) 


00-1 11 




F4(1) 




R3(2) 


1-1100 




Fsd) 




R^(3) 


-1 1000 




F-|(2) 




R2(3) 


0 111 -Cl 




F2(2) 




0 


10 1 1 -C2 




F3(2) 




0 


01 1 1 -C3 




F4(2) 




0 


1 1 0 1 -C4 




F5(2) 




0 


1 1 0 1 -C5 




F-,(3) 




0 






F 2 O) 






F3(3) 






F4(3) 






F5(3) 






^1 






'2 






S3 






S4 






S5 






ALPHA 





Figure A.l. Matrix Representation for Illustration Network. 
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APPENDIX B 



EXPERIMENTAL NETWORKS 





Figure B.2. The 13 Node/60 Link Network. 
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APPENDIX C 

MODEL GENERATION PROGRAM FOR EXAMPLE PROBLEM 

$JO0 1533P 

SQPTIONS FREc 

THIS IS THE DATA GENERATION PROGRAM FOR HPSiil. 

THE FORMAT SHOJLD BE CHANGED ACCORDINGcy. 

N IS NUMBER OF NOCESt L IS NUMBER DF LINKS. 

IEU,L) IS DEFINED AS FOLLOWS: 

1. LINK INCOMING TO NOOEI IS -1. 

2. LINK OUTGOING FROM NODEI IS 1. 

3. LINK UNCONNECTED WITH NOOEI IS 0. 

IZC0L(K,U IS DEFINED AS FOLLOwS: 

l.LINK OUTGOING FROM NODE K IS 0. 

2. OTHERWISE, 1. 

IR(I,KI IS THE INPUT MATRIX FROM NODE I TO NODE K. 
CAP(L,1) IS THE LINK CAPACITY MATRIX. 

SOLNKltJl) IS THE LINK MATRICX CONNECTING NODES. 



DIMENSION INC (13, 60) , IE ( 13 , toO ) , I ZCOL ( 13 , oO) , 

* IRN(13,6J) , IR(13,13J,CAP(60,1) , ICOLNl ( 1 ,60 ) , 

* ICOLNl 13,60), IC0LN2(13, 60) 

INTEGER I ,J,K ,L,N,NJLL, lOEN , KA, J 1 , 1 NC2, 1 A , 

* ID , IQ, IRQWNO,! L ,KC ,KD ,IRN1 

C *** DATA INITIALIZATION ***** 

READ (5, 200 N ,L 

READ (5,225)( I COLNKl ,J1 ),J1 = 1,L) 

READ (5, 22 0) (( IE( I ,J) , 1 = 1 , N) , J=1 , L ) 

READ( 5, 230) ( ( IR( I ,K) , 1= 1 ,N) , K=1 , N) 

READ (5, 2 40) (C AP( I Q,1 ), I Q=l, L ) 

200 FORMAT! 215) 

220 F0RMAT(5I3) 

225 FORMAT! 4110) 

230 FORMAT! 513) 

240 F0RMAT!4E12.5 ) 

C 

C ZERO COLUMN GENERATION **** 

DO 19 J = 1,L 
DO 20 I = 1,N 
INCE = I E! I,J ) 

IFdNCE .EQ. 1) GO TO 21 
INCE = 1 
GO TO 23 
21 INCE = 0 
23 IZCOLd , J) = INCE 
20 CONTINUE 
19 CONTINUE 

C *** COLUMN NUMBER GENERATION 
lAO = 1000000 
DO 1 K = 1,N 
IROWNO = 0 
DO 12 J1 = 1,L 

IC0LN!K,J1) = IC0LN1!1,J1) + K lAO 
IC0LN2!K,Ji) = IC0LN!K,J1J - K 
12 CONTINUE 

C *** INCIDENT MATRICX GENERATION 
DO 2 I = 1,N 
IFd .EQ. K) GO TO 2 
DO 3 J = 1,L 
NULL = IZCOL!K,J) 

IFd .EQ. K) NULL = 0 
lOEN = IE!I,J ) 

INC! I ,J ) = NULL^IDEN 
3 CONTINUE 
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2 

61 

1 

7 

63 

62 

65 

♦ ♦♦ 

66 

67 

♦ #« 

68 

69 

5 

140 

145 

500 

510 

520 

540 

550 

560 

565 

570 

580 



IRQi^NO = I#100 

IRN(I,K) = IRCWNO + K + lAO 

CONTINUE 

ROW SECTION GENERATION 
ID = 1 

DO 61 I = 1,N 

IF(I .EQ. K) GO TO 61 

wRITE(6 ,140)IRN( I ,K) 

CONTINUE 
CONT INUE 

CAPACITY ROW GENERATION 
DO 7 J = 1,L 

WRITE(6,145) IC0LN2(1,J) 

CONTINUE 

COLUMN SECTION GENERATION *** 
WRITE(6,500) 

DO 65 K = i,N 
DO 62 J = 1,L 
NA = 0 

DO 6 3 I = i,N 
IF<I .EQ. Ki GO TO 63 
NULL = IZCJL(K,JI 
IDEN = IE(I,J ) 

INC! I ,J I = NULL* I DEN 
INC2 = INCII,J) 

IF(INC2 .EQ. OJ GO TO 63 
NA = NA + 1 

WRIT E (6, 510)1 COLN (K, J ), IRNI I tK) i INC2 
CONTINUE 

IFINA .EQ. 0) GO TO 62 

WRITE (6 , 520) ICOLNIK, J ) , IC0LN2( K , J ) , I 0 

CONTINUE 

CONT INUE 

CAPACITY GENERATION *** 

DO 6 6 J = 1,L 
lA = J + 1000000 
WRIT£(6,540) lA, ICOL N2( 1 , J) , ID 
CONT INUE 

ALPHA COLUMN GENERATION *** 
WRITE(6,565) ID 
DO 67 J =1,L 

WRIT El 6, 570) IC0LN2( 1 , J ) , CAP I J , 1 ) 

CONT INUE 

RHS GENERATION *** 

WRITE(6,550) 

DO 6 9 K = 1,N 
DO 68 I = 1,N 
IF(I .EQ. K) GO TO 68 
WRITE(6,560) IRN ( I ,.< ) 1 1 R ( I , R ) 

CONTINUE 
CONT INUE 
DO 5 IQ=1,L 
CONTINUE 
WRIT E(6, 580) 

F0RMAT(1X,*E C',I7) 

FORMAT! iX, •£ L* , 17) 

FORMAT! • COLUMNS' ) 

F0RMAT!4X,'X* ,I7,2X, 'C* ,I7,2X,I5 j 
F0RMAT!4X,*X' ,I7,2Xt 'L* fI7,2X,I5) 

FORM AT! 4X,' S' , 17 , 2X, ' L' , 1 7, 2 X, I 5 ) 
FORMAT! ' RHS' ) 

FORMAT! 4X, ' IN PUT ' , 5X , 'C ' , I 7 , 2X, I 5) 
FORM AT! 4X,* ALPHA' ,5X , 'OBJ' , 7X, 15 ) 

FORM AT! 4 X, 'ALPHA' ,5X , ' L ' , 17 , 2X , E 12. 5 ) 
FORMAT! ' ENDAT A' ) 

STOP 

END 
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$ENTRY 

5 a 

010200 020100 01 0300 030200 

020400 020500 040300 050400 



1-10 0 0 
-i 1 0 0 0 

10-100 
0-1100 
010-10 
0 1 0 0-1 

0 0-110 
000-11 
0 0 0 0 0 

0 0 0 0 0 

0 4 0 0 0 

4 0 0 0 0 

0 0 0 0 0 

- 10 . 
- 10 . 

SENTRY 



- 10 . - 10 . 

- 10 . - 10 . 



- 10 . 

- 10 . 
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APPENDIX D 



OPTIMIZATION PROGRAM, INPUT DATA 
MATRIX PICTURE AND OUTPUT FOR EXAJ^IPLE PROBLEM 

^/CHANG JOS ( 1538, 13081 CHANG JOON GONG * ,CL ASS=A 
//* NPSIIi EXAMPLE PROBLEM 

//* 

// EXEC MSSMPS 
//CPC.SYSIN 00 * 

PROGRAM ( • NO* 1 
INITIALZ 

TITLE (‘EXAMPLE PROBLEM, 1ST RUN FOR SUCC-SAT*) 
MOVE IXOBJ,'OBJ*i 
MOVE (XRHS,* INPUT* ) 

MOVE (XOATA, • LINEQS* ) 

MOVE (XPBNAME, 'PROJECT* I 
CONVERT ('SUMMARY*) 

SETUP (*MIN*) 

BCOOUT 

PICTURE 

WHIZARO 

PRIMAL 

SOLUTION 

EXIT 

PEND 



/* 

//EXEC.SYSIN DO « 
NAME LINEQS 

ROQS 
N OBJ 
E C1000201 
E C1000301 
E ClOOOAOl 
£ C1000501 

E C1000102 
E C1000302 
E C1000402 
E Cl 0005 02 
E C1000103 
£ C1000203 

£ Cl 000403 
E C1000503 
E Cl 0001 04 
E Cl 0002 04 
E C1000304 
E Cl 0005 04 
c C1000105 
E C1000205 
E Cl 0003 05 
E C1000405 
E L1010200 
E L1020100 
E L1010300 
E LI 030200 
E L1020400 
E L1020500 
£ LI 0403 00 
E L1050400 
COLUMNS 



X1020101 C1000201 1 

XIO 2 O 1 OI L1020100 1 

X1030201 C1000201 -1 

X1030201 C1000301 1 

XI 030201 L1030200 1 

X1020401 C1J00201 1 

X1020401 ClOOOAOl -1 

X1020401 L1020400 1 
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X1020501 C1000201 

X1020501 C1J00501 

X1020501 L1020500 

XI 040301 C1000301 -1 

XI 0403 01 Cl 000 401 1 

XI 0403 01 L 1040 300 1 

X1050401 C1000401 -1 

X1050401 C1000501 1 

X1050401 L 1050400 1 

X1010202 C1000102 1 

X1010202 L1010200 1 

X1010302 C1000102 1 

XI 010302 CiOOO-302 -1 

X1010302 L1010300 1 

XI 0302 02 Cl 000 302 1 

X1030202 L1030200 1 

X1040302 C1000302 -1 

X1040302 C1000402 1 

X1040302 L1040200 1 

X1050402 C1000402 -1 

X1050402 C1J00502 1 

X1050402 L1050400 1 

X1010203 C1000103 1 

X1010203 C1000203 -1 

X1010203 L1010200 1 

XI 0201 03 C1000103 -1 

X1020103 C1000203 1 

XI 0201 03 LI020100 1 

X1010303 C1000103 1 

X1010303 L1010300 1 

X1020403 C1000203 1 

X1020403 C1000403 -1 

X1020403 L1020400 1 

X1020503 C1000203 1 

X1020503 C1000503 -1 

X1020503 L1020500 1 

X1040303 C1000403 1 

X1040303 L1040300 1 

XI 050403 Cl 000 403 -1 

XI 050403 Cl 000 503 1 

X1050403 L1050400 1 

XI 0102 04 C1000104 1 

X1010204 C1000204 -1 

X1010204 L1010200 1 

XI 0201 04 C 1000 104 -1 

X1020104 C1000204 1 

X1020104 L1020100 1 

X1010304 C1000104 1 

X1010304 C1000304 -1 

X1010304 L1010300 1 

XI 0302 04 C 1000204 -1 

X1030204 C1000304 1 

X1030204 L1030200 1 

X1020404 C1000204 1 

X1020404 L1020400 1 

X1020504 C1000204 1 

X1020504 C1000504 -1 

XI 02 05 04 L1020 500 1 

XI 05 0 4 0 4 C1000504 1 

X1050404 L1050400 1 

X1010205 C1000105 1 

X1010205 C1000205 -1 

X1010205 L1010200 1 

X1020105 C1000105 -1 

X1020105 C1Q00205 1 

X1020105 L1020100 1 

X1010305 C1000105 1 
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X1010305 


C1000 305 


-1 




X1010305 


LiOlOSOO 


1 




X1030205 


Cl 000 205 


-1 




X1030205 


C1000205 


I 




X1030205 


L 1030 200 


1 




X1020405 


CiJ00205 


1 




X1020405 


C1000405 


-1 




X1020405 


L1020400 


1 




X1020505 


C1000205 


1 




XI 0205 05 


L 1020 500 


1 




XI 040305 


Cl 000 305 


-1 




X1040305 


C1000405 


1 




X1040305 


LI 040 300 


1 




SlOOOOOl 


L1010 200 


1 

1 




S1000002 


L1020100 




S1000003 


L1010300 


1 




SI 0000 04 


L1030200 


1 




SI 000005 


L1020 AOO 


1 




S1000006 


L1020500 


1 




Si 000007 


L1040300 


1 




S1000008 


L1050A00 


1 




AL PHA 


oeo 


1 




ALPHA 


L1010200 


-0. lOOOOE 


02 


ALPHA 


L1020 100 


-O.IOOOOE 


02 


ALPHA 


L1010300 


-O.IOOOOE 


02 


ALPHA 


L1030200 


-0. lOOOOE 


02 


AL PHA 


L1020400 


-O.IOOOOE 


02 


ALPHA 


L1020500 


-O.IOOOOE 


02 


AL PHA 


L1040 300 


-0. lOOOOE 


02 


ALPHA 

RHS 


L1050400 


-0 .lOOOOE 


02 


INPUT 


C1000201 


0 




INPUT 


C1000 301 


0 




INPUT 


C1000401 


0 




INPUT 


C1000501 


0 




INPUT 


Cl 000 102 


0 




INPUT 


01000 302 


Q 




I NPUT 


C1000402 


0 




INPUT 


C1000 502 


0 




INPUT 


C1000103 


0 




INPUT 


C1000203 


5 




INPUT 


C 1000 403 


0 




INPUT 


C1000 503 


0 




INPUT 


C1000104 


2 




INPUT 


C1000204 


0 




INPUT 


C1000 304 


0 




INPUT 


C1000504 


0 




INPUT 


C1000105 


0 




INPUT 


C1000 205 


0 




INPUT 


CIOOO 305 


0 




INPUT 

ENOATA 


C1000405 


0 





// 
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tXAMPLE fROaLtM, 1ST fUN FOR SuCCfcSSIVt SAT 
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tAAMPLfc fHOBLEM, 1 ST fUN fCK SUCCtSSlVE SAT 
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EXAMPLE FfiOflLEMi i^C FUN FCR SUCCESSIVE 5AI 
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EXAMPLE FP 06 LEH, iNC fUN FOR SUCCESSIVE SAT 
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EXAMPLE FRGaLtM, FUN fCR SUCCESSIVE SAT 
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EXAMPLt PROBLEM, 2 PC PUN FOR SUCCESSIVE SAI 
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IXAMFLE fROBLEM, ^ IH FUN fCH SUCCfcSSlVt SAT 
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tXAMPLE PBCaLEM, <,lh fUN fCR SUCCtiSIVE SAT 
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EXAMPLE EMCdLEM, SIP PUN fCR SuCCESilVc SAT 
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EXAMPLE EPOBLEM, tih EON fCR SUCCESSIVE SAI 
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tXAMPLE fHCULErt, tlh H.N KiH iUCCLiiiVE SAT 
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IXAMPLfc FhOBLEM, tU FUN FDR SUCCfcSSiVt SAI 
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APPENDIX E 



OPTIMIZATION PROGRAM, INPUT DATA 
I4ATRIX PICTURE AND OUTPUT FOR MAX-SLACK APPROACH 



//CHANG J38 ( 1538, 1808) CHANG JOON wOONG • ,CLASS = A 

//* 

//# MPSIII =XAWPL£ PfiCBLcM 

//« 

// EXEC MSSMPS 
//CPC.SYSIN DD » 

PROGRAM (• NO* ) 

INITIALZ 

TITLE {'EXAMPLE PROBLEM, 2ND RUN FOR MAX-SLACK*) 
MOVE <XOBJ,*OBJ*) 

MOVE {XRHS,* INPUT*) 

MOVE (XDATA, • LIMEQS* ) 

MOVE (XPBNAME, 'PROJECT* ) 

CONVERT ( * SUMMARY * ) 

SETUP {'MAX') 

8CD0UT 

PICTURE 

WHIZARC 

PRIMAL 

SOLUTION 

EXIT 

PEND 

/* 

//EXEC.5YSIN DO ♦ 

NAME LINEQS 

ROWS 



N OBJ 
E Cl 0002 01 
E C1000301 
£ C1000401 

E Cl 0005 01 
E C1000102 
£ C1000302 

E C1000402 
E C1000502 
E C1000103 
E Cl 0002 03 
£ C1000403 

E C1000503 
E C1000104 
E Cl 0002 04 
E C1000304 
£ C1000504 

E C1000105 
E C1000205 
E C1000305 
E C1000405 
E L1010200 
E L1020100 
£ L1010300 

E L1030200 
£ L1020400 

E L1020500 
E LI 040300 
E L1050400 
COLUMNS 



X1020101 C1000201 1 

X1020101 L1020100 1 

X1030201 C1000201 -1 

X1030201 C1000301 1 

X1030201 L1030200 1 

X1020401 C1000201 1 

XI 020401 C 1000 401 -1 
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XI 020401 L1020400 1 

X1020501 C1000201 i 

X1020501 C1000501 -1 

X1020501 L1020500 1 

X1040301 C1000201 -1 

X1040301 ClOOO^Ol 1 

XI 0^03 01 LI 040 300 1 

X1050401 C1000401 -1 

X1050401 C1000501 1 

X1050401 L1050400 1 

X1010202 C1000102 1 

X1010202 L1010200 1 

X1010302 C1000102 1 

X1010302 C1000202 -1 

X1010302 L1010300 1 

X1030202 C1000202 1 

X1030202 L1030200 1 

XI 040302 C1000 302 -1 

X1040302 C1000402 1 

XI 0403 02 L1040 300 1 

XI 050402 C 1000402 -1 

X1050402 C1000502 1 

X1050402 L1050400 1 

X1010203 C1000I03 1 

X1010203 C1000203 -1 

X1010203 L1010200 1 

X1020103 C1000 103 -1 

XI 0201 03 C1000 203 1 

X1020103 L1020100 1 

X1010303 C1000103 1 

X1010303 L1010300 1 

X1020403 C1000203 1 

X1020403 C1000403 -1 

X1020403 L1020400 1 

XI 0205 03 C 1000 203 1 

X1020503 C1000503 -1 

X1020503 L1020500 1 

XI 0403 03 C 1000 403 1 

XI 040303 L1040200 1 

X1050403 C1000403 -1 

XI 050403 C1000503 1 

XI 050403 L1050400 1 

X1010204 C1000104 1 

X1010204 C1000204 -1 

X1010204 L1010200 1 

X1020104 C1000104 -1 

X1020104 C1000204 1 

X1020104 L1020100 1 

X1010304 C1000104 1 

X1010304 C1000204 -1 

X1010304 L1010200 1 

X1030204 C1000204 -1 

X1030204 C1000204 1 

X1030204 L1030200 1 

XI 020404 C1000 204 1 

X1020404 L1020400 1 

X1020504 C1000204 1 

X1020504 C1000504 -1 

X1020504 L1020500 1 

X1050404 C1000504 1 

X1050404 L1050400 1 

X1010205 C1000105 1 

X1010205 C1000205 -1 

X1010205 L1010 200 1 

X1020105 C1000105 -1 

X1020105 C1000205 1 

X1020105 L1020 100 1 
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* 

* 

* 

* 

* 

RHS 



X10103G5 
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INPUT 
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I NPUT 
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INPUT 
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2.5 







ENDATA 

/* 

// 



1 

1 

1 
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APPENDIX F 







LINK FLOW 


FOR 9/36 NETWORK 








Table 


F.l . Link 


Flow for 9/36 Network 


( ) ; MAX-SLACK 


LINK 




DESTINATION NODE 




IOTA! 












1 U 1 r\L 


(NODES) 


1 2 


3 4 


5 6 7 


8 


9 FLOW 


1-2 


0.5 
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3.5 




(4) 

4 
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1-3 
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4 
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1-9 
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4-3 


1 .75 








1 .75 


3-5 






(4) 

3.75 




(4) 

3.75 
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4-6 
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LINK 




DESTINATION NODE 








TDTAl 


(NODES) 1 
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3 4 5 6 
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FLOW 
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Wl 
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/ sJ 
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REQ-MAT 




Rl(5)=15, R5(2)=9 
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LINK FLOW FOR SYMMETRIC REQ-MAT OF 13/60 NETWORK 
Table H.1. Link Flow for Symmetric Req-Mat of 13/60 Network ( ) ;MAX-SLACK 
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REQ-MAT Ri(13)=2, Ro(12)=2, R.(9)=2, Rq(5)=2, R,^(2)=2, RiJ1)=2 
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